Gauge theory description of glass transition 
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An analytical approach, which develops the gauge model of the glass transition phenomenon, is 
suggested. It is based on the quantum field theory and critical dynamics methods. The suggested 
mechanism of glass transition is based on the interaction of the local magnetization field with 
the massive gauge field, which describes frustration-induced plastic deformation. The example of 
the three-dimensional Heisenberg model with trapped disorder is considered. It is shown that the 
glass transition appears when the fluctuations scale reaches the frustrations scale, and the mass of 
the gauge field becomes equal to zero. The Vogel-Fulcher-Tammann relation for the glass transition 
kinetics and critical exponent for non-linear susceptibility, 1.7 < 7 < 3, are derived in the framework 
of the suggested approach. 

PACS numbers: 



The formulation of a universal theory of glass tran- 
sition has remained one of the most intriguing but still 
unresolved problems of condensed matter physics 
for a long time. Many systems which manifest this phe- 
nomenon regardless of their nature enable us to conclude 
that this phenomenon does not depend on any micro- 
scopic details, but is determined by the symmetry prop- 
erties of the systems, as in the case of phase transitions, 
regardless of the scale. However, it is well known that 
the glass transition is not accompanied by divergence of 
the static correlation length of the local magnetization 
and linear susceptibility, \l , typical for critical phenom- 
ena. This does not allow to assign this glass transition 
phenomenon to phase transition. To support the opinion 
that glass transition is a critical phenomenon, it should 
be noted that critical slowing down of all processes takes 
place at the microscopic level during glass transition and 
the nonlinear susceptibility, xn, diverges. Besides, re- 
cent computer simulations enable us to suggest that glass 
transition is still a specific form of critical phenomena, 
in which the scale of a temporary ordered region of the 
medium-range crystalline order is considered to be a di- 
verging static correlation length ||. These are the fea- 
tures of a critical phenomenon, however, the reason for 
this still remains unclear. 

The theoretical approach to the glass transition de- 
scription, introduced in this paper, is based on well 
known methods of the critical phenomena description 
and nonequilibrium dynamics, and promotes the ideas 
mentioned in ^-|TT[. According to this theory the glass 
transition is accompanied by the formation of local order- 
ing in the same way as the second order phase transition. 
We suppose that the key feature, which distinguishes 
one phenomenon from the other, is the presence of the 
disorder- induced frustrations. The frustrations can be 
induced by both quenched randomness, associated with 
the random sign of the spin interaction in the case of 
the spin-glass transition, and topological pro per ties of 
the system, as the geometrical frustration in []10|. The 
frustrations distort the system's space B , make it hetero- 
geneous, and lead to its plastic deformation. We describe 



these frustrations below as topological defects (discli- 
nations), which are sources of appropriate gauge fields 
f§ §, The fluctuations of the local magnetization 

field in the disordered phase increase as the tempera- 
ture approaches the phase transition point, T c . However, 
the system contains the frustrations, whereas the ordered 
domains should not contain them. Therefore, the fluctu- 
ations grow until their size remains less than the typical 
spacing between the frustrations, Lfi uct Lf rust . Since 
the fluctuations can not grow more than L f rus t , the criti- 
cal slowing-down occurs while they reach this size. Below 
we show that this critical behavior is attributed to van- 
ishing the effective mass of the gauge field, which hap- 
pens at the temperature T g > T c . As a result, the sys- 
tem freezes in a state with a disordered structure. Thus, 
the glass transition represents the critical slowing-down 
of fluctuations growth by disorder- induced frustrations, 
that complies with the "frustration-limited domain the- 
ory" (FLDT) jnj and correlates with other approaches 

In more detail, let us consider an example of the three- 
dimensional static isotropic Heisenberg model, whose or- 
der parameter, s, is considered to be a local magnetiza- 
tion vector with the SO (3) continuous symmetry group. 
The pure model describes the system undergoing the 
paramagnetic-ferromagnetic phase transition, with no 
quenched disorder available. The static action of the 
model is well known and has the form of 
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-(9,s) 2 + C/(s) 



dr. 
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where U(s) = /i 2 s 2 +?;s 4 , . Note that at /i 2 < the action 
is invariant under the SO (3) gauge transformations and 
(s) = 0. However, at /i 2 < the symmetry is explicitly 
broken, since the system can arbitrarily "choose" only 
one state from all equivalent states with the minimum en- 
ergy U(s) potential minima situated on the |s| = ifi/y/2v 
sphere. We fix the vacuum by means of fixing a point 
on the sphere. The system is no longer symmetrical with 
respect to the SO (3) gauge group, but it is invariant un- 
der the SO(2) group of the rotation around the chosen 
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direction. 

In a disordered spin system the equilibrium spin orien- 
tations at different points are not collinear and a suitable 
connection between orientations is done by introducing 
a gauge field, A", and replacing the ordinary derivative, 
diS, by a covariant derivative, DiS, [0, [l2|: 



S = 



where 
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dr, (2) 



(3) 
(4) 



J" is the source of the Af. field which is introduced to 
the action in addition to the general part. The reason 
for this is discussed below. 

The method of introduction of disorder in the theory 
is very important and plays a key part. In order that 
a spin system possess glass properties we should inject 
to this system a disorder which gives rise to frustration 
of its structure. Unfortunately, before ones did not pay 
attention to this fact properly. In ||, 0, for example, 
the gauge field is free, and any quenched randomness is 
absent. This is right when the dynamic soliton solutions 
appear spontaneously in the form of pair of oppositely 
directed vortexes. However, in our case the frustrations 
should be quenched and static. In || and M the disorder 
is associated directly with the quenched gauge field, and 
the A a ^ field is frozen in an arbitrary configuration with 
some P(A") distribution function. It is supposed that the 
quenched gauge field describes the frustrations. But it is 
not quite right too, since presence of the quenched gauge 
field in the system does not mean yet presence of frustra- 
tions. The quenched frustrations can be represented by 
immovable disclinations passing through the frustration 
planes ||, |ll) (see Fig. 1). They are the static sources of 
the A® gauge field, and should be injected in the model 
by means of the static sources field, J°. Therefore, in 
contrast to |j, we believe that it is more correct to 
consider the source field, J", but not the gauge field, 
in the capacity of quenched random field. In this case 
the A° field remains to be dynamic one. For illustra- 
tion the A® field can be interpreted as a local relative 
rotation of neighboring spins, which corresponds to their 
local equilibrium. Spins can be movable, but their local 
equilibrium configuration around a static frustration in 
the SV volume, bounded by the SS sphere, should satisfy 
to the following condition: 



\§F%dS v = J J^O, 

SS SV 
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which follows from the principle of least action. 

Expansion of the local magnetization field, s, near one 
of the vacuum states, for instance (s) = (0, 0, i[i/V2v), 
in small <fi = s — i[i/\/2v deviations, and use of the 



gauge transformation properties allow to rewrite (£]) in 
the form of the action of two massive vector bosons, Af, 

(k = {1, 2}), with the mass Mq = ig/i/V2v, one massless 
vector boson, A^, and one scalar field, </>: 
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Note, that this Lagrangian is not gauge-invariant in pres- 
ence of external sources, that is important below. 

Let us consider the system with randomly quenched 
frustrations, which are described by means of the sources 
field, J°. The J° field is static and absolutely random. 
Therefore, for simplicity, we assume that ( J°(r) J°(r')) = 
IoS(r — r'), where Iq is the intensity of the quenched dis- 
order. It is shown below that this parameter is propor- 
tional to some "structural temperature" . Averaging over 
leads to redefinition of the partition function: 
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exp 
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Dct>DAl, 
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where J . . . Dx is the continual integral. It leads to ad- 
ditional contribution to the Af. "mass" , which takes the 
form of 



(8) 



Thus, the frustrations lead to the renormalization of the 
gauge field mass. 

The renormalization of the gauge field mass affects 
the critical behavior of the system, since it shifts 
M = singularity to the temperature range above the 
paramagnetic-ferromagnetic transition point, T c . If we 
assume that fi 2 — aks{T — T c ), where a is some con- 
stant, then from (||) we have the critical divergence of the 
A^ field correlation radius at T g = T c + AIqv / 'aksg 2 in 
the paramagnetic phase. One can suppose that this can 
lead to critical slowing-downs of the fluctuations. Thus, 
the disorder-induced frustrations inhibit the growth of 
the cf> field correlation length, and the system freezes in a 
disordered state, that conforms with the glass transition 
description in FLDT. 

In order to examine the above assumption we should 
investigate the kinetics and susceptibility of the consid- 
ered model near T g . To do this, we will consider the 
non-equilibrium dynamics of the system near the criti- 
cal point. This is a classical problem. Therefore, one 
can apply either the method of the dynamic generating 
functional , or the classical limit of the Keldysh tech- 
nique JlJ] , since both these methods coincide in the clas- 
sical limit. We will choose units such that ksTg = 1. 
The usage of the functional technique for the description 
of non-equilibrium dynamics leads to the representation 
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of the partition function of the system in the form of 



where 

S* = 
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dvdt, (10) 



4> = {<p, </>}, and A^ = {A 1 ^, are vectors, the com- 
ponents of which are named "quantum" and "classical" 
respectively [Q, G _1 and A^J are matrices, inverse to 
the Green functions matrices having the following form: 
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Note that opposed to Yang-Mills theory here there is no 
necessity to introduce an additional ghost field in the 
theory because of the gauge symmetry breaking and the 
presence of the gauge field mass. All Green functions of 
the theory are explicitly determined by the same reason. 
The components of the Green function of the scalar field 



G* A \k,u) = 
G K (k, u) 
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where is the kinetic coefficient, which corresponds to 
the local magnetization. The components of the Green 
function of the massive gauge field are: 
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where Ta is the kinetic coefficient, which corresponds to 
the gauge field. Their graphic forms are presented in 
Fig. 2 

The interaction of the local magnetization fluctuations 
per gauge field plays the key part in the considered 
theory (Fig. 2). This becomes clear in the process of 
renormalization. Close to M 2 — the gauge field be- 
comes massless, but the scalar field remains massive with 
/i 2 = AI^v/g 2 . Therefore, the contribution to the renor- 
malization is made only by the loops of the gauge field 
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FIG. 1: a) is the retard part of A^u field Green function , 
A/Li b) is the Keldysh part of this Green function, A^,; c) 
is the retard part of the Green function of the scalar field <f>, 



G ; d) is the Keldysh part of this function, G 
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FIG. 2: 



propagators. According to the separation of massive field 
theorem (TH) the Feynman diagrams, containing the prop- 
agators of field, the mass of which is appreciably larger 
than the external momentum, are inversely proportional 
to the degree of this mass, and make a finite contribution 
to the renormalization. In Fig. 3 some graphs giving log- 
arithmically divergent contributions to the renormalized 
theory are presented. It is easy to check that near T g the 
considered theory is renormalizable, since the logarith- 
mically diverging corrections lead to the renormalization 
of only existing terms of the Lagrangian. However, it 
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FIG. 3: 

is necessary to consider in detail the contribution to the 
renormalization of the node of the a diagram, which 
is shown in Fig. 3. This term is interesting because the 
node is proportional to the relaxation time of the lo- 
cal magnetization field, and determines the kinetics of 
the glass-transition: the loop of the light field, A® which 
is given in Fig. 2, makes the logarithmically divergent 
contribution ~ ln(l/A)5(w), where A is the regulariza- 
tion parameter (cutoff of momentum) . This term deter- 
mines the divergence contribution of the a graph (Fig. 3) 
in which the massive field loop leads only to multiplying 
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the logarithm by Al^v/g 2 factor: 



^MVA) 1 e -^\t\/r^ dt = 



(16) 



where £ Q is the time of the observation of the system. 
One can see that for Jo — > 0, or with a short observation 
time, T^g 2 /AIqv 3> £ , this contribution becomes negli- 
gibly small, and the theory becomes nonrenormalizable. 
This relates to the symmetry properties of the Yang-Mills 
model. In this case the fluctuation dissipation theorem 
(FDT) is always broken because of free energy trans- 
fer between the local magnetization modes and massless 
gauge field modes (Goldstone modes). However, there 
is also some problem in the presence of the quenched 
disorder: this node becomes nonlocal in time, which vi- 
olates the correctness of the renormalization procedure. 
In order to avoid this problem it is possible to divide 
the time space of the system into two intervals: when 
t *C r c/,g 2 1 AIqv this contribution is negligibly small; but 
when t a 3> T ^g 2 / AI$v this graph makes the logarithmi- 
cally divergent contribution to the renormalization. 
Hence in the one-loop approximation the renormalization 
group has the form: 
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where £ = ln(l/A), e — 4 — d. Actually, this form of the 
renormalization group takes into account the FDT viola- 
tion on small time scales, and seems to be natural for the 
quasiergodic system. In |]l6| it was shown, that match- 
ing the marginal solutions in ( |l7j ) leads to the Vogel- 
Fulcher-Tammann relation for the temperature depen- 
dence of the system relaxation time near T g . From the 
condition of the stable point existence, d\n(g 2 )/d(, — 0, 
d\n(v)/d£ — 0, we get g 2 = ir 2 s, v — g 2 /2, and, using 
the result of lEfll, 



Trel = OC exp 
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Hence, the critical slowing down of all relaxation pro- 
cesses does occur in T g , and follows the Vogel-Fulcher- 
Tammann relation. 

It is evident from (|l2|) that the linear susceptibility, 
XL = d{4>)/dh ~ \i~ 2 = g 2 /AIov (h is an external source 



of the field 



is finite in T„ 



yj g 2 /AIqv, is finite too. 



The correlation length, 
However, nonlinear sus- 
ceptibility, xn = d 3 (<fi)/dh 3 [|l7| , diverges near T g : it is 
not difficult to check that the nonlinear susceptibility (see 
Fig. 4) encloses infinitely increasing contributions. The 
simplest one of them corresponds to the h diagram in 
Fig. 3, and gives the divergent contribution proportional 
to (e£)~ 3! /( 2+3e / 2 '. Hence in one loop approximation it 



is possible to estimate xn oc (T — T g ) 7 for T 
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where 12/7 < 7 < 3. The lower limit is determined by 
the critical dimension, d = 4, the upper limit by the real 
dimension, d — 3. This is in good agreement with the ex- 
perimental observations |17 . Thus, one can assert that 



the system freezes in the state with a disordered struc- 
ture of the local magnetization field, and T g is the glass 
transition temperature. 
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FIG. 4: Graphic representation of the nonlinear susceptibility, 
Xjv = (4>4><j>(j>) fc=o (the external momenta are k = in the 
momentum representation) 

In conclusion it is necessary to note once again, that 
the gauge symmetry breakdown is not spontaneous in the 
considered model. It is induced by the random sources 
of the gauge field, which is connected with the topologi- 
cal frustrations. The symmetry breakdown leads to the 
appearance of the gauge field mass. However, the gauge 
field mass becomes small when the temperature reaches 
T g (T g > T c ). Then the scale of the gauge field fluctua- 
tions becomes commensurable with the scale of the inho- 
mogeneities induced by the frustrations. Close to T g the 
local magnetization fluctuations are small, but the gauge 
field fluctuations dramatically increase, which leads to 
the critical slowing-down of the gauge field dynamics, 
as well as the dynamics of the local magnetization field. 
This slowing-down freezes the disordered structure of the 
local magnetization field with the finite correlation length 
and susceptibility. If frustrations are absent, Iq — ► 0, the 
freezing temperature coincides with the phase transition 
temperature, and the scales of the ordered domains be- 
come macroscopic. Then the diagrams with the loops of 
<j> and 4> fields become divergent, and the system experi- 
ences the paramagnetic-ferromagnetic phase transition, 
which is described within the standard critical dynam- 
ics @. 

It should also be noted, that it is not only quenched 
randomness that leads to the appearance of the gauge 
field mass. The disorder in the system can also be in- 
duced by the topological frustrations. The geometrical 
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frustrations connected with the topological features of 
the condensed matter structure can produce a massive 
gap of A't. For example, in the model of the defected 
states of the orientation order || the mass of the gauge 
field is determined by the summed topological charge, Q, 
which is connected with the curvature of the hypersphere, 



corresponding to the locally preferable structure. 
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